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ABSTRACT 

Context. We study an SDSS sample of galaxies within ~50 Mpc of the nominal center of the Ursa Major supercluster. 

Aims. Our aim is to study galaxy distribution around groups in the supercluster and the link between the distribution of relaxed and 

nonrelaxed galaxy systems with respect to the supercluster environment. 

Methods. Using the FoF algorithm, 40 galaxy groups were identified in this region. Velocity distributions for these groups were 
studied after applying a 3cr-clipping routine for outlier removal. We classified the systems according to the results of normality and 
substructure tests applied to member galaxies. Then, we studied the relative distribution of relaxed and nonrelaxed systems across the 
supercluster. 

Results. We find that 68% of galaxy groups are Gaussian and that all the non-Gaussian systems have substructures and probably 
correspond to multimodal systems in redshift space. We also find that the Gaussian systems inhabit the denser regions of the super- 
cluster, with higher densities of both red and blue galaxies within 2.5 /r'Mpc, and have smaller group-group pairwise separations. 
Conclusions. Our results suggest a spatial segregation of dynamical states, where relaxed systems may have formed and evolved 
earlier and faster around high-density peaks, while nonrelaxed systems may be growing slower on the peripheries of lower density 
peaks. In this picture, galaxy clustering seems to be prompting a continuous internal evolution in the supercluster, with several groups 
collapsing into the more evolved and contracted regions. 
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1. Introduction 

Groups of galaxies contain most of galaxies in the Universe 
and are the link between individual galaxies and large-scale 
structures (e.g., Huchra & Geller 1982; Geller & Huchra 1983; 
Nolthenius & White 1987; Ramella et al. 1989). The dissi- 
pationless evolution of these systems is dominated by grav- 
ity. Interactions over a relaxation time tend to distribute the 
velocities of the galaxy members into a Gaussian distribution 
(e.g. Bird & Beers 1993). Although the theoretical line-of-sight 
(los) velocity distribution expected in galaxy systems is not ex- 
actly Gaussian (e.g. Merritt 1987; Kazantzidis et al. 2004), phe- 
nomenological evidence has been suggesting for a long time that 
normality can be assumed for systems in dynamical equilibrium 
(e.g. Yahil & Vidal 1977). 

Several works have shown that the analysis of velocity dis- 
tributions can separate galaxy groups into relaxed or nonre- 
laxed systems (e.g. Hou et al. 2009; Ribeiro et al. 2010, 2011). 
Actually, several important properties of groups can be stud- 
ied from this perspective. For instance, Hou et al. (2009,2012) 
find rising velocity dispersion profiles for non-Gaussian (NG) 
groups. They have also shown that the majority of NG groups 
have substructures. Ribeiro et al. (2010) found that galaxies 
are redder and brighter in Gaussian (G) groups out to 4R2oo- 

Send offprint requests to: A.L.B. Ribeiro 



Martinez & Zandivarez (201 1) have studied the luminosity func- 
tions of G and NG groups, concluding that G groups have a 
brighter characteristic magnitude. Einasto et al. (2012a, 2012b) 
studied velocity distributions to find signatures of multimodality 
in SDSS clusters. 



The NG nature of the velocity distribution of NG groups 
must be caused by recent group-group mergers. At the same 
time, NG groups could be contaminated by chance projections. 
It is therefore important to establish the real nature of NG groups 
and how they are linked to the environment. We study an SDSS 
(Sloan Digital Sky Survey) sample of galaxies within ~50 Mpc 
of the nominal center of the Ursa Major Supercluster aiming 
to establish a relation between large-scale environment and the 
evolution of galaxy groups. After identifying groups in the su- 
percluster using the FoF (Friends-of-Friends) algorithm, we split 
up the sample into relaxed and nonrelaxed systems according to 
the velocity distribution of member galaxies. Then, we investi- 
gate the distribution of galaxy groups in the supercluster with re- 
spect to their dynamical state. The paper is organized as follows. 
Section 2 presents the data and methods we have used, Section 3 
presents our analysis, and Section 4 briefly summarizes the main 
conclusions of the work. 
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2. Data 

In this work, we have used galaxies from the SDSS-DR6 
(Adelman-McCarthy et al. 2008) within -50 Mpc around the 
coordinates of the cluster AbeU 1383 ( =177.46°, Dec=54.62°, 
and z=0.059) taken here as the nominal (and arbitrary) "center" 
of the Ursa Major supercluster (hereafter UMaS). The redshifts 
are in the range 0.045< z <0.075. We have found 18,204 galax- 
ies in this region with apparent magnitude limit m r = 17.77. 
Clusters were identified using the FoF algorithm (Huchra & 
Geller 1982), following the choice of linking lengths by Berlind 
et al. (2006) [b ± = 0.14 and b\\ = 0.75Q which were optimized 
for SDSS. The FoF algorithm with these linking lengths finds 
galaxy groups with N ga \ >10 that have an unbiased multiplicity 
function (see Berlind et al. 2006). These settings led us to find 
40 galaxy groups in the UMaS, with numbers of galaxies in the 
range 15 < A^ ga i < 94. Our catalog contains all but one of the 23 
galaxy systems known in the region[^](see Kopylova & Kopylov 
2009). 

Velocity distributions for the identified groups were stud- 
ied after applying a 3<x clipping routine for removing outliers. 
First, a virial analysis [^]is done following Carlberg et al. (1997), 

where fl 2 oo = V3~cr/[10i/(z)] and M 20 o = 3R 20 oo- 2 /G. Then, 
we applied three well known normality tests to the cleaned ve- 
locity distributions - the Anderson-Darling (AD), Shapiro-Wilk 
(SW) and Jarque-Bera (JB) tests - to decide if a galaxy sys- 
tem is Gaussian or NG. Our criterium is that at least two tests 
should reject normality at the 95% confidence level to be con- 
sidered NG; if not, it will be considered Gaussian. Finally, we 
applied two substructure tests - the /3 (West et al. 1988) and A 
tests (Dressier & Schectman 1988) - to decide if a galaxy system 
has substructures. Our criterion is that a galaxy system is consid- 
ered to have substructures if at least one of the two tests rejects 
regularity with 95% confidence level. Results are presented in 
Table 1, where columns are (1) group identification; (2) and (3) 
equatorial coordinates in degrees; (4) average redshift of mem- 
bers after removing outliers; (5) velocity dispersion in km/s; (6) 
R 2 oo radius in /r'Mpc; (7) logarithm of M 2 qo in h~ l M Q ; (8) total 
number of galaxies after outlier removal; (9) results for the NG 
nature of the los velocity distribution: for Gaussian and 1 for 
NG distributions; (10) results for substructure tests: for regular 
and 1 for irregular systems; (11) Abell cluster identification. 
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Table 1. Properties of 40 galaxy groups identified in the UMaS. 
Columns 9 and 10 indicate results for normality and substructure 
tests described in the text. 



3. Analysis 

3.1. Multimodality 

We find that 27 out of 40, i.e. 68%, of galaxy groups are 
Gaussian and that all the NG systems have substructures. Also, 
we find that four groups have substructures (as indicated by the 
P test) but have Gaussian velocity distributions according to our 
criteria. If we include them as nonrelaxed systems, the total frac- 
tion of virialized groups would fall to 58%. We take a conserva- 
tive view to define nonrelaxed groups and assume that substruc- 
tures are dynamically significant only if they are also associated 
to non-Gaussianities in the velocity distribution of a galaxy sys- 
tem. 

Following Ribeiro et al. (201 1), we applied the Kolmogorov- 
Smirnov (KS) and the Cramer-von Mises (CvM) two-sample 



1 The parameters b ± and b\\ are the projected and line-of-sight linking 
lengths in units of the mean intergalaxy separation. 

2 The missing cluster is beyond our search radius. 

3 In this work, cosmology is defined by £l m = 0.3, C1 A = 0.7, and 
H = lOO/i km s-'Mpc-'. 



tests on the G and NG subsamples. We find that NG velocity 
groups are signifcantly larger and more massive at the 95% con- 
fidence level for both tests, with best fit relations: 



R™ a (1.27 ± 0.22) R° 



21 Ml 



- M NG 



(2.00 ± 0.65) M% 



200 



with [(p=0.0065 - KS), (p=0.0019 - CvM)] for the radius distri- 
butions, and [(p=0.0029 - KS) (p=0.0017 - CvM)] for the mass 
distributions. This result is similar to the one found by Ribeiro et 
al. (2011), a possible consequence of applying the virial analy- 
sis to nonvirialized objects. Actually, we can assume the velocity 
distributions of NG groups as Gaussian mixtures with unknown 
number of components, and use the Dirichlet process mixture 
(DPM) model to study the velocity distributions (see Ribeiro et 
al. 201 1 for details). The number of components in the mixture is 
computed using the R language and environment (R Develoment 
Core Team) under the dpmixsim library (da Silva 2009). We find 
the following results: 

- Unimodal groups: 3, 9, 34 
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Fig. 1. Four examples of multimodality diagnostics for groups in the UMaS region. Vertical bars depict individual galaxies in each 
group. The solid lines indicate the best Gaussian mixture for for radial velocities distributions in km/s. The dashed lines show the 
individual components of the multimodal systems. The dotted lines indicate the average Gaussian curve for 1,000 resamplings of 
groups 03 and 34, unimodal but non-Gaussian systems. 



- Bimodal groups: 4, 6, 18, 19, 22, 25, 26, 28, 30 

- Trimodal group: 1. 

In Figure 1 we plot the velocity distributions for typical cases 
of NG groups. Multimodality seems to be the prominent cause 
of non-Gaussianity in galaxy group velocity distributions, corre- 
sponding to 10 out of 13, i.e. 77% of the NG cases in the UMaS. 
There are three unimodal systems (#03, #09, and #34) all hav- 
ing negative kurtosis excess (k = -1.26, -0.87, and k = -1.05, 
respectively) - the mean values for Gaussian and NG groups are 
kc - -0.34 and &ng = —1.16 A negative kurtosis corresponds 
to a flattened and usually multimodal distribution. However, if 



the separation between the modes is too small, the DPM model 
is unable to recognize distinct components. In Fig. 1, we il- 
lustrate two cases of NG unimodal systems, also plotting the 
Gaussian curve resulting from 1,000 resamplings keeping 90% 
of the dataj^For each group we apply the normality tests men- 
tioned in Section 2 to the new velocity distribution. The reduced 
samples classified as normal distributions are used to compute 
the average Gaussian we plot in Figure 1 . 



4 For instance, group #03 has 26 members, thus each resampling pro- 
vides a reduced group with 23 members. 
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Fig. 2. Distribution of galaxies and galaxy systems in the UMaS in supergalactic coordinates. Galaxies are in black; Gaussian groups 
and non-Gaussian are in red and blue, respectively. 



3.2. Supercluster Environment 

Now, we try to link the dynamical state of galaxy systems and 
the supercluster environment. The spatial distribution of galaxy 
groups in the UMaS region is presented in Figure 2 in su- 
pergalactic Cartesian coordinate^ where we see galaxies and 
groups. From this figure, it seems that G groups inhabit the 
denser regions of the supercluster, approximately following the 
central arch-shaped structure in the plot. That could be hinting at 



5 The supergalactic coordinates are based on the preferred plane of 
the Local Supercluster. After the coordinate transformation, the Abell 
cluster A1383 was chosen to be at the center of galaxy distribution. 



some environmental influence in the galaxy groups spatial dis- 
tribution. 

We have used two simple indicators to probe the possible re- 
lation between galaxy groups and the supercluster environment. 

3.2.1. Distribution of galaxies around groups 

First, we built the stacked G and NG galaxy number density (in 
redshift space) up to 2.5 /z _1 Mpc, around the groups' centers, a 
distance corresponding to 6/?2oo for our sample. Results are 
presented in Figure 3, where we see similar shapes for the pro- 
files around G and NG groups, but with densities systematically 
higher in the case of G groups. Error bars in this plot are obtained 
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Fig. 3. Galaxy number density profiles in redshift space. Filled 
and open circles denote G and NG systems, respectively. Error 
bars are obtained from 1,000 bootstrap resamplings of galaxy 
groups. Top, middle and bottom plots refer to all, red, and blue 
galaxies, respectively. 



from a bootstrap technique with 1,000 resamplings applied to 
the galaxy groups. This result indicates that G systems present 
significantly richer surroundings than NG sytems, with ratio of 
number densities of galaxies around G and NG groups following 
log[vNG(/)/vG(/)] — 0.35 ±0.10 (vis the galaxy number density). 

We also classified galaxies into blue or red objects following 
the criterion of Strateva et al. (2001) - a red galaxy should have 
u - r > 2.22. This enables us to rebuild the galaxy density pro- 
files for these two populations. The results for red galaxies are 
also presented in Figure 3 (the middle plot), where we see that 
G groups have more red galaxies within 2.5 /r'Mpc than NG 
groups log[v^ G (r)/v 1 G (r)] 0.46 + 0.11, but note that the dif- 
ference are smaller at larger radii, reflecting (i) the well known 
correlation between colors and environment - the fraction of red 
galaxies decreases toward the outer part of clusters (e.g. Balogh 
et al. 2000; De Propis et al. 2004), and (ii) the quicker decline 
for G systems. At the same time, for blue galaxies the difference 
between G and NG systems is smaller, but still statistically sig- 
nificant: log[v^ G (r)/v G (r)] a; 0.69 ± 0.12, see Figure 3 (on the 
bottom plot). These results show that G groups are denser and 
redder than NG groups up to 2.5 Zi 'Mpc. 



3.2.2. Distribution of groups through the UMaS 

Also, we studied the relative separation between groups in both 
G and NG subsets. We computed the distance from each group to 
the first (k = 1), second (k = 2), and third (k = 3) nearest groups 
in the spatial distribution. In Figure 4, we illustrate our results for 
the supergalactic YZ projectiorj^] The mean pairwise distances 
of G groups show just a small variation with k: 11 /i^'Mpc (k = 
1), 12 /r'Mpc (k = 2), and 14 h^Mpc (k = 3); while we see a 
wider variation for NG groups: 10 /z _1 Mpc (k = 1), 19 /i^'Mpc 
(k = 2), and 26 /!~'Mpc (k = 3). This suggests that G groups are 
closer to each other than NG groups. 

We estimated the significance of the pairwise separation test 
by scrambling the identities of the G and NG groups 10,000 
times, preserving the overall numbers of 27 Gaussian and 13 
non-Gaussian. We calculated the number of replicas with more 
extreme statistics than the input data (i.e. larger differences be- 
tween G and NG mean pairwise separations). For k = 1, we 
found 4,356 cases, indicating that G and NG systems have the 
same distributions with respect to the first neighbor. However, 
for k - 2 and k — 3 we found 279 and 232 cases of more ex- 
treme statistics. We can therefore reject the null hypothesis of the 
same spatial distribution for G and NG systems at the ~2.8% and 
~ 2.3% signficance levels, respectively. The difference between 
the spatial distributions of G and NG groups is reinforced by 
the result of the two-dimensional KS test (e.g. Press et al. 1992; 
Fasano & Franceschini 1987), which evaluates the null hypoth- 
esis that a pair of data sets are drawn from the same distribution. 
We found that this is rejected at the 95% confidence level for all 
2D data projections. 

Combined, these results indicate that G groups probably in- 
habit denser regions, with higher densities of red galaxies, and 
are closer to each other (at least for k - 2 and k — 3) than NG 
groups - a trend that suggests an important environmental effect 
on structure formation through the UMaS volume. 



We found approximately the same results for the XY and XZ pro- 
jections 
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Fig. 4. Pairwise distances from each group to the first (k — 1), second (k = 2) and third (k = 3) nearest neighbors in the supergalactic 
YZ projection. The arrows point to the kth nearest neighbor among groups of the same class. The red and blue arrows indicate 
distances for G and NG systems, respectively. The mean distance between groups appears on the right side of each box. The size of 
all boxes is 100 x 100 h~ 'Mpc. 



4. Discussion 

Superclusters are the largest structures identified in the universe. 
They correspond to complex systems comprising a few up to 
several galaxy clusters and groups experiencing long-term gravi- 
tational evolution. In this work, we studied a galaxy sample ~ 50 
Mpc around an arbitrary center (the Abell cluster 1383) of the 
Ursa Majoris Supercluster. We have addressed two questions: 
galaxy distribution around groups in the supercluster and the link 
between the distribution of relaxed and nonrelaxed galaxy sys- 
tems with respect to the large-scale environment in which they 
are embedded. We find that 68% of galaxy groups are Gaussian 



and that all the NG systems have substructures. This corresponds 
to a lower fraction of relaxed groups in comparison to our pre- 
vious studies (Ribeiro et al. 2010, 201 if] This result cannot 
be due to the different criteria employed to split up G and NG 
groups, since we are being more conservative to define nonre- 
laxed systems in the present work. Alternatively, this may re- 
sult from the supercluster environment, which could be favoring 
interactions and mergers between the groups. Actually, finding 
multimodality as the main cause of non-Gaussinity seems con- 



But consistent with Hou et al. (2009) for CNOC2 data. 
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sistent with the idea of ongoing processes in nonrelaxed systems 
(e.g. Ribeiro et al. 2011; Hou et al. 2012; Einasto et al. 2012a). 

In this respect, an important issue relates to the possibility 
of also having NG groups as the result of spurious members in 
the velocity distributions. Wojtak et al. (2007) show that most 
statistical methods (including the 3<x-clipping procedure used in 
the present work) exclude 60%-70% of unbound galaxies on av- 
erage. The probability of chance alignments with the remaining 
interlopers was estimated with the following procedure: 

1. We assume that all groups are 100% reliable in the member- 
ship. 

2. We place all galaxies within them at the center, making their 
velocity dispersions nil. 

3. And we rotate the cube around a random axis. 

4. Then, we disperse the los velocities according to a Gaussian 
distribution for each group with the original velocity disper- 
sion for the group. 

5. Finally, we apply an FoF code with the same parameters 
mentioned in Section 2, and analyze the Gaussianity of the 
resulting groups. 

After 1,000 repetitions of this procedure, we found that 
~15% of the mock groups, which are Gaussian by construction, 
become NG, once additional members are included. This means 
that about two NG groups in our sample would correspond to 
Gaussian systems with spurious members in their velocity dis- 
tributions. This fraction is too low to signicantly modify the re- 
sults we found in the previous sections, so most of NG groups 
are probably portraits of multimodal systems. 

Even presuming this is so, Gaussian groups are still domi- 
nant across the supercluster, which suggests the system is not 
young enough. Actually, Einasto et al. (2012b) show that su- 
perclusters have two main morphological types, spiders and fil- 
aments. Clusters in superclusters of spider morphology have 
higher probabilities to have substructure and higher peculiar ve- 
locities of their main galaxies than clusters in superclusters of 
filament morphology. The Ursa Major supercluster is classified 
as filament-type in that work, which suggests it might be dynam- 
ically evolved. In the same direction, we found that G groups in- 
habit denser regions, with higher densities of both red and blue 
galaxies within 2.5 /i~'Mpc and are closer to each other than 
NG groups, suggesting that relaxed systems may have formed 
and evolved earlier around high-density peaks, while nonrelaxed 
systems may be growing more slowly in the peripheries of lower 
density peaks. This shows the complex nature of superclusters 
so that further studies are needed to understand the role of the 
large-scale environment in galaxy group distribution and evolu- 
tion. 
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